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Final Report 

The Prediction of Noise Scattered by a Blended Wing-Body/Ducted Fan Configuration 


Introduction 

In this proof of concept research, a computational method was developed for predicting the 
sound field created by the scattering of ducted fan engine noise by a blended wing-body (BWB). It was 
assumed that all acoustic processes were linear and time harmonic with excitation frequency to Inflow 
effects were neglected and no penetration boundary conditions were applied to the engine nacelle and 
BWB surfaces. 

A scattering approach was adopted in which the total acoustic field is written as the sum of 
known incident (from the engine duct) and unknown scattered parts. We further assume that the 
incident field is independent of the scattered field Application of the above conditions to the equations 
of linearized acoustics yields the Helmholtz equation (reduced wave equation) for the scattered pressure 
with Neumann boundary conditions. 

Research Approach and Results 

An equivalent source method (ESM) was used to solve the classical Helmholtz boundary value 
problem. The ESM was chosen for its ease of implementation and computational simplicity relative to 
the more accurate boundary integral or boundary element methods. See references 1-3 and attachments 
1-2 of this document for ESM details. The ESM solution methodology development proceeded in four 
stages. 

1) 2-D Scattering: The ESM was applied to scattering of 2-D, point source generated incident sound by 
simple 2-D elliptical objects. Numerical studies were designed for this simple problem so that pertinent 
discretization parameters could be identified and applied to the 3-D scattering problem. 

2) Simple 3-D Scattering: The 2-D ESM was extended to include the 3-D scattering of point source 
generated sound by an ellipsoid. Attachment 1 contains sample results and a description of the ESM 
implementation. All calculations were performed on a Pentium II 350MHz processor with 256Mb 
RAM. Numerical experimentation showed that meaningful results were obtained for values of the non- 
dimensional parameter a < 250 Where a = kW, k = oj/c, c = ambient sound speed, W = largest spatial 
dimension of scattering body We estimate that the maximum value of a can be doubled by employing 
the latest PC technology (as of 1 1/99) and doubled again with advanced algorithm development 

3) Scattering of Ducted Fan Noise by BWB geometry: Using drawings of the BWB design, the BWB 
planform shape was approximated by deforming and connecting several ellipsoids. This approach yields 
an analytical representation of the scatterer that is simple to evaluate. Point source generated sound for 
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the ESM was replaced by incident sound from a simple ducted fan model. The ducted fan noise 
prediction program TBIEM3D (reference 4-5) was used to calculate the incident sound Sample 
acoustic field calculations for a 4-bladed fan were obtained Details of the implementation and 
computational results appear in attachment 2. 

4) Comparison with Langley BWB Shielding Experiment: Due to discretization size limitations 
discussed in stage 2 above, a direct comparison between the ESM and the experiment is impossible We 
show that shielding results obtained with the ESM agree qualitatively with the BWB measurements To 
simplify calculations we considered a smaller excitation frequency, larger duct radius, and shorter duct 
length than the experiment. Incident noise was generated by a point monopole placed in the center of 
the middle duct. Radiation from the outer two ducts was ignored for this set of calculations. By 
considering a shorter center duct and neglecting the two outer ducts, TBIEM3D computational time is 
much reduced. These factors do not influence the ESM discretization and could be included at the cost 
of computational time. See table 1 below for a description of the pertinent operating parameters 



ESM 

Experiment^ 

kR 

6.0 

11.5 

W/R 

36.7 

69.4 

Ld/R 

2.0 

9.74 

LAV 

0.49 

0.51 

HAV 

0.14 

0.084 

a = kW 

232 

798 


Table 1 : Fundamental Parameters 
R = duct radius L D = center duct length W = BWB wingspan 
L = BWB length (tail to nose) H = BWB maximum thickness 
(*) Experimental data obtained from Clark and Gerhold 

In accordance with the BWB experiment, sound pressure levels were calculated on a plane 
parallel to the BWB planform and located 0.29 wingspans beneath the duct axis. Computations are 
performed with and without the BWB presence. For contrast we also compute SPL’s on a plane 0 29 
wingspans above the BWB The ESM results are shown in figure 1 and clearly show the scattering 
effects due to the BWB 

In figure 2 , the difference in sound pressure levels between the shielding and no shielding cases 
are plotted for both the upper and lower observation planes. The results in figure 2a show a shadow 
region beneath the BWB that is somewhat similar in shape and magnitude to that produced by the 
experiment. 
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It is noted that the ESM was not designed to take advantage of the symmetry present in the BWB 
experiment. Doing so would approximately double the maximum attainable a. Therefore, applying 
symmetry arguments and the latest PC technology to the ESM would nearly quadruple the maximum 
a meeting the conditions of the BWB experiment for 1XBPF excitation frequency For higher 
harmonics, ESM computational requirements exceed PC capabilities. 

Conclusions 

Due to its simplicity and ease of implementation, the ESM was shown to be an effective noise 
prediction tool for acoustic scattering problems. Only simple scattering bodies were considered, but the 
method is applicable to arbitrary shapes. Comparison with BWB experimental results showed that the 
ESM calculations are physically correct. 

In essence, the ESM is an approximation to boundary integral or boundary element methods. 
The computational advantages of boundary methods over finite difference and finite element methods 
are significant and have been addressed in reference 5. Like any solution technique requiring 
discretization, the ESM is limited to moderate values of excitation frequency and body size. 

To expand the range of applicability of the ESM, the consumption of computational resources 
must be reduced. Theoretical innovation is required for this purpose. For example, the use of multipole 
expansions combined with boundary integral and advanced grid generation techniques might be used to 
lower the number of equivalent sources required to produce the desired boundary conditions 
Furthermore, ESM algorithms are well suited for multiprocessor computing which would result in a 
significant decrease in computational time. 
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SPL, dB 



34 36 38 40 42 44 46 48 50 52 54 



Figure 1 : ESM Acoustic Field Calculated on Observation Planes 
0.29 Wingspans Above and Below Duct Axis k = 6.0 m' 1 
a) no shielding b) BWB present 



Figure 2: Calculated Noise Difference Due to BWB k = 6.0 m’ 1 
a) lower plane b) upper plane 
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X, inches 


Calculated overall engine-radiated noise 
reduction by the Blended Wing Body in flyover. 


Figure 3: BWB Experimental Results - Obtained from C. Gerhold 
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Attachment 1 : Report on the Development of a Noise Prediction Method for a Blended Wing/Ducted 

Fan Configuration 



Report on the Development of a Noise Prediction 
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Scattered by a Blended Wing/Ducted Fan Configuration”. 





Incident sound generated 
by a time harmonic, monopole 
point source is scattered by 
an ellipsoid 






Equivalent Source Method 
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ESM is much simpler to implement and calculate. 




Placed at Mesh Vertices 
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least squares methods. 



Results 



source and collocation points. 















Re [Instantaneous Acoustic Pressure] 







Numerical experimentation with this problem has provided 
guidance for conducting the final phases of this research project. 





Estimated Completion Date: 09/30/1999. 











More Future, Future Effort 
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Attachment 2: Application of the Equivalent Source Method to the Scattering of Ducted Fan Engine 

Noise by a Blended Wing-Body 
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Outline 





Incident sound generated 
by a time harmonic noise 
source is scattered by 
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mose mat use boundary integrals or boundary elements, the 
ESM is much simpler to implement and calculate. 
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Comparisons with experimental results are pending. 



Oblique View 









D Scattering Results - Axial Planes of Sound Pressure Level 






-D Scattering Results - Horizontal Planes of Sound Pressure Level 





Conclusions 



ESM algorithms are well suited for multi-processor computing. 




